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Abstract 

The recently published data for for J ftp production in Pb-Pb collisions at 158 A 
GeV are analyzed. For low values of transverse energy Et, where normal suppression 
dominates, (p1)(Et) scales with the path length of the gluons which fuse to make the 
J ftp. In the Et domain of anomalous suppression (pl)(Et) is found to rise linearly with 
the relative amount of anomalous suppression. This empirical law is reproduced within 
an analytically solvable transport model which allows high pt J ftp's to escape anomalous 
suppression. Interpreted in this way, the data for (p1)(E t ) lead to an estimate of ~ 4 
fm/c for the duration of anomalous suppression. 

New data mean new surprises. This has been a recurrent phenomenon in the study of 
charmonium suppression in high energy nuclear collisions during the last decade of research. It 
also holds for the recently released data of transverse momentum distributions for charmonia 
taken in Pb-Pb collisions at 158 A GeV 0. In particular for central collisions where anomalously 
large suppression has been observed, the new data display interesting features: the mean values 
(pi) of transverse momentum p t increase slowly with transverse energy E t and turn steeply 
upward where the J ftp yield drops (see Fig. 1). Can this behavior tell us something new about 
the nature of anomalous J ftp suppression? 

The previous data for the transverse momentum dependence of J ftp production in nuclear 
pA and AB collisions and their interpretation are presented in two reviews [0, ||. Basically, 
two mechanisms are proposed: 

(i) Rescattering of gluons in the initial state |§, ||, |]: In a pA collision, the gluon of the 
projectile proton scatters from target nucleons before it fuses with a gluon from the 
target to form the J ftp. Gluon rescattering in the initial state is treated as a random 
walk in transverse momentum and the observed (pi) is predicted to increase linearly with 
the mean length £ g of the path of the incident gluon. In an AB collision both gluons which 
fuse to the J ftp are affected by the rescattering. This effect in the initial state has been 
clearly identified in the data, also for AB collisions HJ. For central AB collisions it had 
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Figure 1: Data for the mean transverse momentum {pf) of J/^'s produced in Pb — Pb collisions 
at 158 A GeV as a function of transverse energy E t together with the data for the relative J/i/j 
production cross section S = a(Pb + Pb — > J/ip + X)/a(Pb + Pb — > DY + X) in the same 
experiment [0]. The range of E t values is divided into the domains of normal suppression ("N") 
(E t < 50 GeV) and anomalous suppression ("A") (E t > 50 GeV). Within the domain 'A", the 
subdomain "fl" with E t > 100 GeV indicates events of high E t which rise from fluctuations. 
The solid line SN(E t ) is calculated for normal nuclear suppression "N"[[U|, while S Q b denotes 
the observed values. 
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been predicted |7], [8j that initial state rescattering together with anomalously strong Jf tp 
absorption in the final state may lead to a saturation or even a decrease of £ g and therefore 
of (pi) for the very large values of E t . The new data for Pb-Pb collisions contradict this 
prediction. 



(ii) Escape of high p t J ftp's in the final state P, ITTDL ITTI: Within the scenario of the quark-gluon 



plasma (QGP) anomalous charmonium suppression occurs in a limited space-time region 
during final state interaction. Only, J ftp's with sufficiently high transverse momenta 
Pt have a chance to escape the "deadly" region. Therefore anomalous suppression acts 
preferentially on low pt J ftp's and the surviving charmonia should show higher values of 
(pi) with increasing amount of anomalous suppression. This is what indeed is seen in the 
new data. Do the data then confirm this mechanism? 

In this letter we analyze the new data for (pi) of J ftp production in Pb-Pb collisions and 
identify two empirical laws: (a) In the domain of normal J ftp suppression the values for (pi) 
depend linearly on the path lengths of the gluons. (b) A linear correlation between the values 
of (pi) and the relative amount of anomalous suppression is discovered in the E t range where 
anomalous suppression has been identified. To our knowledge the second empirical law has 
not been noticed before. We discuss various explanations and then present an analytically 
solvable transport model based on mechanism (ii), i.e. escape of high p t J ftp's from the region 
of anomalous suppression. A more detailed analysis of all the data from [0J in the light of both 
mechanisms (i) and (ii) is in preparation. 

Fig. 1 shows the data of ref. [|1|] for (p1)(E t ) of J ftp's as a function of transverse energy 
E t together with the previously published data on J ftp suppression in the form of the ratio 
S(E t ) = a(Pb + Pb^ J ftp + X) fa(Pb + Pb -> DY + X). We divide the E t region into three 
domains: (a) Small transverse energy (E t < 50 GeV), where there is no anomalous suppression 
and the yield is well described by the Glauber approach with an effective Jf tp nucleon absorption 
cross section (which correctly describes pA and S-U collisions). We denote quantities in this 
domain by the index "JV" (for "normal"), (b) For values E t > 50 GeV one observes anomalous 
suppression, i.e. the data S Q b(E t ) deviate from the the predictions of the Glauber approach. 
We denote all quantities in this E t region by the index Ll A" (for "anomalous"), (c) Within the 
anomalous region the data at very high values, E t > 100 GeV, show a particular behavior: the 
suppression drops while the data for (pi) rise. These high values of E t correspond to the most 
central collisions and are interpreted [|T2| , p~3[] to arise from fluctuations in the transverse energy. 
We have indicated this origin by the symbol at the appropriate places. 

We begin our analysis of the data by investigating to which degree (for which values of E t ) 
the mechanism (i), gluon rescattering in the initial state, explains the data. This mechanism 
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leads to a dependence 



^) AB (Et) = (p 2 t ) NN + ^if(E t ), (1) 



where (p 2 ) N is the contribution already present in J/ip production in the elementary NN 
event, while the second term is linear in the mean length £ g (E t ), which the two gluons travel 
in nuclear matter before they fuse. The constant in front of £ g depends on (p 2 ) gN , the mean 
transverse momentum acquired in a gluon-nucleon collision, and \ 9 n, the mean free path of 
a gluon in nuclear matter. This constant is taken as an adjustable parameter, while £ 9 {E t ) is 
calculated. Then the observed values of (p 2 ) AB (E t ) are plotted versus the calculated values 
of £g(E t ). If a linear relation appears, one takes it as support for the gluon rescattering 
mechanism. This analysis has already been done in ref. [p]] for the Pb-Pb data. However, 
£ AB (Et) is calculated neglecting absorption of the J/ip in the final state. We therefore repeat 
their analysis with absorption. Of course, the results depend on the employed absorption model. 
We choose the one of ref. |13 , which correctly reproduces the data for J/ip suppression even 



in the domain of fluctuations. Then 



/ d 2 b d 2 s dz A dz B (if(s, z A ) + £ B (b- s, z B )) K(b, s, z A , z B , E t ) 



iAB {Et) = - '— ; v \\ ir j ~ /y „v : - - - , (2) 

where 



9 J d?b d 2 s dz A dz B K(b, s, z A} zb, E t ) 



£ A (s,z A ) = / dzp A (s,z)/p , 
J ~ 00 

£ B (b-s,z B )= \ dzp B (b - s,z)/p (3) 



'-D 



are the geometric lengths which the two gluons travel through the nuclear density along the 
^-direction with impact parameter b, and p(r) is the density of nuclear matter. To simplify the 
numerical calculations, we use in the following the uniform distribution. The expression for the 
kernel K in eq. @ is taken from ref. fl~3|] : 



K(b,s,z A ,z B ,E t ) = p A (s,z A )p B (b- s,z B ) 




2D 



dzp A (s,z)+ / dzp B (b — s,z) 



r n p (bJ)\-P(E t \b). (4) 
\Et){b) ) 

Here, the final state absorption for the J/ip is contained in the exponential ("normal" absorption 
by nucleons with an absorption cross section cr^[f) and in the Theta function ("anomalous" 
suppression suddenly setting in, when the density of produced matter exceeds a threshold 
density n c ). The function P(E t \b) describes the distribution of transverse energy in events at a 
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Figure 2: The observed values (pf) for J/ t/>'s from Pb—Pb collisions plotted versus £ g (E t ) , 
the mean path of the gluons before they fuse to the J/ip. The values i g (E t ) are calculated within 
the model |TB[ which contains normal and anomalous nuclear suppression and which is able to 
reproduce the cross section data in Fig.l. As defined in Fig.l we have divided the E t range 
into domains "N", "A" and "fl": Only the five data points on the straight line belong to the 
domain of normal suppression, while all data on the backward branch belong to anomalous 
suppression. 



given impact parameter b. We have followed ref. [13| in the notation and the numerical values 
for the constants and therefore refer the reader who is interested in more details, to this paper. 

Fig. 2 shows a plot of the experimental values of (pf)(E t ) versus the calculated values of 
£^ B (E t ). The qualitative picture is rather different from the one given in ref. [I]: Only those 
values of E t , which correspond to the domain of normal ("/V") suppression in Fig. 1, show a 
linear relation. We fit it by a straight line and find a value for the slope constant 



l gN 



{pi) 



gN 



\gN 



(0.102 ± 0.006) {GeV/cf/fm, (5) 



which differs from the value (0.081±0.003) (GeV/c) 2 /fm given in ref. since we use a different 
prescription for the calculation of ig. In the E t domain where anomalous suppression ("A") 
sets in, the plot shows an anticorrelation: With increasing values of E t , the calculated values 
of tgiEt) decrease, while the experimental values for (pf)(E t ) increase. This behavior in the 
anomalous domain hides some physics other than gluon rescattering and will be studied in the 
following. 

We assume that all J/^'s which are suppressed by an anomalous mechanism have still 
another source influencing the transverse momentum distribution above the one from gluon 
rescatering in the initial state. If we decompose the observed cross section S b(E t ) for J/ip- 
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production (relative to .DF-production) into its two contributions 



S ob (E t ) = S N (E t )-S A (E t ), (6) 

where the contribution Siy{Et) from normal absorption is a theoretical quantity and is calculated 
by using the kernel K from eq. ([J), without the Theta function and is shown in Fig. 1 by the 
solid line. The anomalous suppression, SA(E t ), is defined as the difference between observed, 
S b(E t ), and calculated values, SN(E t ). We associate different values, (p^)jv(-E't) and (p 2 )A{E t ), 
with the normal and anomalous contributions, respectively. Then the observed value (p 2 ) b(E t ) 
can be written as 

S N (E t )(p 2 ) N (E t ) - S A (E t )(p 2 t } A (E t ) 



(pt)o b (E t ) 



S b(E, 



(Pt)N(E t ) + ^^M(E t ), (7) 



where 

5p 2 (E t ) = (p 2 t ) N (E t ) - (pl) A {E t ). (8) 

Eq. (0) should be valid in the full E t domain comprising anomalous suppression and normal 
one (where Sa{Ei) = 0, by definition). For normal suppression, the E t dependence of the 
observed values (p 2 ) b(E t ) is solely carried by values (p^)jv(-E't) calculated from eqs. (fj-fl) 
without the 0-function and displayed in Fig. 2. The calculation of (p 2 )N{E t ) can also be 
extended into the domain of anomalous suppression (E t > 50 GeV). We find that {pj)N{E t ) 
changes by maximally 3% over the whole domain 0.25 < Sa/S i, < 1.25, where data exist and 
are displayed in Fig. 3. Therefore, in the domain of anomalous suppression, the E t dependence 
of {p 2 ) b(E t ) is predominantly carried by the second term in eq. (0). In order to disentangle 
the E t dependence residing in Sa/So^ from the one in Sp 2 (E t ), we plot the observed values of 
(p 2 )ob(E t ) versus the ratio SA(E t )/ S b(E t ). The result is displayed in Fig. 3. To a very good 
accuracy all points lie on a straight line with a slope 

dp 2 = (0.132 ± 0.007)(GeV/c) 2 . (9) 

The result indicates that the E t dependence of (p 2 ) b(E t ) in the domain of anomalous suppres- 
sion arises predominantly from the ratio SA{E t )/ S b(E t ). We also draw the attention to the 
E t region where fluctuations dominate: While a sudden drop in the J/ip suppression and 
a sudden rise in (p 2 ) occur (Fig. 1), the empirical regularity shown in Fig. 3 continues into 
the domain without any noticeable change in character indicating that no new physics 
appears in this domain of E t . The empirical correlation displayed in Fig. 3 is the first new 
result of this paper. 
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Figure 3: The observed values (pf)(E t ) (in the "A" domain, E t > 50 GeV) for J/-0's from 
Pb—Pb collisions plotted against the relative amount of anomalous suppression S^Et)/ S Q b(E t ). 
Here, SA(E t ) is the difference between the calculated value Sjy(E t ) and the observed values 
S a b(E t ) in Fig.l. (The data in the "N" domain, not shown, would all be found at Sa/ 'S Q b = 0). 

Can we understand this relation? Are we able to calculate the constant 5pf7 We will discuss 
two possible explanations, intermediate x production and escape of high p t charmonia in the 
final state and start with the x- 

The observed intensity of J/ip's produced in nuclear collisions has two contributions, directly 
produced J/^'s and those arising from produced ip h s and x's which decay into J/ip long after 
the collision but before detection. The indirect contributions amount to about 40% |fj, §]. Since 
ip n s and x' s are l ess bound than the J /if), it has been argued || that anomalous suppression 
should act predominantly on the x's and ip n s rather than on the directly produced J/^'s. The 
x's may have a smaller value of (pf) x than the J/ip and ip', since the x' s can be produced 
directly by fusion of two gluons. J / ip and ip' production involves three gluons, where the third 
one is presumably radiated off from a color octet intermediate state and generates additional 
(pf). If anomalous suppression acts predominantly on the Xi the linear relation eq. (0) seems 
plausible with 

% 2 ~ {(pD n jH - (pX n ) ■ ( 10 ) 

Unfortunately, we have no calculations for the difference in eq. (|10"D . Furthermore, the emission 
of a gamma (x — > J/ip + 7) reduces the difference by about 0.1(GeV/c) 2 and may even change 
the sign of bp\. If the mechanism of intermediate x production was the dominant explanation 
for the observed value of bp\ , the values of (pf) observed in the production of f/>"s should have no 
anomalous values. However, the contrary is true according to the data in Pb-Pb collisions Jl|. 
Therefore the contribution of x to the observed J/ijj does not seem a compelling explanation 
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for the behavour of (p 2 ) in the anomalous region. 

We come to the second explanation: Only high p t charmonia escape anomalous suppression. 
This argument is not new|| [10|, [TT[| . In this paper we start from this idea and propose a solvable 
transport model, in order to see whether the empirical relation Fig. 3 and the value 5p 2 , eq. 
(H), can be understood. 

We introduce the phase space distribution n(f,p,t) for J/ip's produced in a Pb-Pb collision 
with a definite value of E t . In a system where the J/^'s have zero longitudinal momentum, the 
values of r*and p denote the transverse position and momentum of the J/ip's, respectively. We 
denote by t = the time, when all collisions involving nucleons (production N + N — > J/ip + X 
and suppression J/ip + N — > no J/ip) have ceased. Then at t = we have what we call normal 
suppression and therefore have 



S N = J d 2 fd 2 pn(f,p,0) , (11) 
(p 2 t ) N = [ d 2 fd 2 pp 2 n(f,p,0)/S N (12) 



for the intensity SV and (p 2 )n of the J/ip, respectively. Here and in the following we consider 
one particular impact parameter (one value of E t ). For t > the phase space distribution 
n(f,p,t) evolves under the influence of anomalous suppression. In order to obtain analytical 
results, we have to introduce three simplifying assumptions: (a) Instead of the absorption 
process happening continuously for t > until all charmonia are "eaten up" , we let anomalous 
suppression happen at one particular time t = t^, while for < t < t^, the phase space 
distribution evolves freely, i.e. 

n(f,p,t) = n((f+ vt),p,0) , (13) 

where v = p/M is the velocity of a J/ip with momentum p within a non-relativistic approx- 
imation, (b) Anomalous suppression happens in such a way that all J/ip's, whose position 
r is smaller than a given radius are suppressed (where depends on E t , rA{E t )). This 
assumption corresponds to the G function, eq. ([|), introduced in ref. [0. Then 



S ob (E t )= d 2 rd 2 pn(r,p,t A ) (14) 

Jr<r A (E t ) 



(p 2 t ) ob (E t )= d 2 rd 2 pp 2 n(r,p } t A )/S ob (E t ) (15) 

Jr<r A (E t ) 

for the observed intensity S ob (E t ) and the observed (p 2 ) h(E t ), respectively. With the further 
assumption (c) of a Gaussian for the initial phase space distribution 

(r 2 p 2 \ 



8 



all integrals in eqs. (|IT| - p^) can be evaluated explicitly and after some algebra one arrives at 



{ptU(E t ) = (p 2 t ) N + In (l + • 5p 2 , (17) 

s P 2 t = (vDn ■ ^ftlR- 2 - (18) 

Eq. (|17D for the mean squared transverse momentum (p 2 ) b(E t ) in the domain of anomalous 
suppression is the second new result of our paper. It shows the dependence of this quantity on 
the degree of anomalous suppression S A = Sjy — S a b, and one recovers the empirical law, eq. (|7|) 
after expanding the logarithm for small values of S A /S b (the logarithmic dependence seems to 
be an artifact of the Gaussian shape for the phase space distribution eq. (jnf)). Furthermore, 
eq. (|T^) gives an analytical expression for the value 5p 2 . It involves the mean squared transverse 
velocity of the produced J/ip, the mean squared transverse radius R of the overlap zone of the 
two colliding nuclei (for central collisions, R is related to the mean squared radius of the Pb 
nucleus) and the time t A , when anomalous suppression happens. Using the empirical values for 
dp 2 = 0.132(GeV/c) 2 , eq. (|), (p 2 t ) N = 1.77(GeV/c) 2 and R 2 = \ (r 2 ) Pb , one finds t A = 4 fm/c. 
This value is not unreasonable in view of values of 5 — 7 fm/c discussed in the literature [|i~4|j 
for the lifetime of the fireball in central Pb-Pb collisions. We expect our value for t A to change 
somewhat, if some of the above assumptions (a-c) are relaxed. However, the basic dependence 
of 5p 2 on the quantities R, (p 2 )n and t A should remain (also for dimensional reasons). 

In summary, we have been able to show, that normal nuclear suppression goes along with 
values of (pf) which depend linearly on the path lengths £^ B of the gluons which fuse to form 
the J/ip, Fig. 2. In the domain of anomalous suppression, the data show a linear dependence on 
the relative amount of anomalous suppression S A /S a b. This relation also holds for the domain of 
fluctuations in E t and supports the interpretation that no qualitatively new mechanism begins 
for these very high values of E t . The empirical relation in the domain of anomalous suppression 
is consistent within the geometric picture "high ptJ/ip's escape anomalous suppression", which 
has been proposed || [H| and worked out JllJ long ago. 
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